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Structural properties of even-even nuclei along the -stability line are systematically analyzed and 
investigated by means of our improved potential-energy-surface (PES) approach. The equilibrium 
deformations obtained from the PES minima are compared with available experimental data and 
other theoretical results. From a brand-new angle of view, taking the nuclei on the stability line 
as footholds, we systematically present two-nucleon separation energies and half lives in the corre- 
sponding isotopic chains. Very regular laws can be seen, especially below Z = 84. The observed 
trends and anomalous behaviors are briefly discussed. The energy difference between theory and 
experiment for the binding energies and two-nucleon separation energies indicates that the improved 
PES calculations with our fitting macroscopic model parameters are relatively in better agreement 
with experiments. Based on the systematic law of two-nucleon separation energies, the bound limit 
of even-even nuclei is crudely estimated along the stability line. 


PACS numbers: 21.10.Re, 21.60.Cs, 21.60.Ev 


I. INTRODUCTION 


Most of the nuclei in nature were formed at the cre- 
ation of the universe or after supernovae explosions [1]. 
Their half-lives are usually longer than the expected age 
of the Earth and the Solar System (~ 4.6 billion years) 
or the Universe (~ 13.8 billion years), indicating they are 
stable or practically stable. These stable nuclei, less than 
300 nuclides (the nuclide number is slightly different, de- 
pending on time scale used), form a ‘valley of stability’ 
(which is also referred as ‘belt of stability’ or ‘band of 
stability’, or a ‘peninsula of stability’ in a ‘sea of insta- 
bility’) in the nuclear chart. The nuclei in the valley are 
of special interest in nuclear physics because they are not 
only naturally existent but also the important candidates 
of the projectile-target combinations in nuclear reactions 
used for producing radioactive nuclei though the high- 
intensity radioactive beam facilities and the radioactive 
targets start to be used nowadays [2]. 

Along the valley of stability, the semi-empirical ap- 
proximation formulas were fitted half a century ago, 
such as N — Z = 0.4A?/(A + 200) in Ref. [3] and 
Z = A/(1.98 + 0.0155A?/3) in Ref. [4], which can allow 
one to estimate the mass number, A, or neutron number, 
N, for a given proton number, Z. The curves described 
by such two formulas are usually referred to as the 8- 
stability lines (also, the backbone of the belt of stability) 
and almost overlap except for the slightly different ex- 
trapolations in the high-Z region. The nuclei located at 
the northwest (southeast) of the stability line contain too 
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many protons (neutrons) to be stable, tending to become 
stable ones by, e.g., a-, G~-, 8*-decay or electron capture 
and allowing them to move closer to the optimal neutron- 
proton ratio. As is well known, only certain number of 
combinations of neutrons and protons can form the stable 
nucleus which has enough ability (e.g., binding energy) to 
permanently hold its nucleons together (however, regard- 
less of neutron-proton combinations, no stable members 
are identified in some isotopes with Z > 83). It has been 
found that the nuclei consisting of specific nucleon num- 
bers, e.g., 2, 8, 20, 28, 50, 82 and 126 (for neutrons), have 
prominent stability, possessing a relatively higher aver- 
age binding energy per nucleon. Such magic numbers 
corresponding to shell closures have been predicted by 
the nuclear shell-structure model [5-9] and proved by ex- 
perimental observables with sudden discontinuities [10], 
e.g., in the nucleon separation energies with changing Z 
or N. 


Up to now, the heaviest doubly magic nucleus is 
298Pbi26, almost lying on the stability line. Naturally, 
the search for next heavier doubly magic nucleus and 
the northeast limit of the nuclear chart is more fasci- 
nating. Theoretical models based on the macroscopic- 
microscopic approaches [11-20] using various potentials 
predicted that the heavier doubly magic candicate is 
the superheavy nucleus ??§Flig4, which lies very closely 
to this extrapolated 6-stability line and may be stable 
enough to exist in nature. Some self-consistent mean- 
field calculations prefer Z = 120 and N = 172 to be the 
shell closures [20-23]. Furthermore, around the next dou- 
bly magic nucleus, the predictions also show that there 
may exist a batch of shell-stabilized superheavy nuclei 
which form a small so-called “island of stability” (relative 
to the “peninsula of stability”). It is expected that nuclei 


near the top of the island of stability can exist for, e.g., 
many years instead of some milliseconds [24]. Presently, 
along the stability line, the experimentally synthesized 
heaviest even-even nucleus is the nucleus 785Noi62 [25] 
though the periodic table of elements has contained 118 
known and named elements [26-29]. 


As mentioned above, nuclear structural properties will 
strongly affect the stability of nuclei. Nuclear binding 
energies and decay half-lives increase markedly in the 
vicinity of nucleon closed shells which give rise to spher- 
ical nuclear shapes. Indeed, nuclear shape is very help- 
ful for understanding the fundamental interaction related 
to microscopic nuclear structure. It is known that nu- 
clear stability both in the superheavy and drip-line re- 
gions depends on nuclear deformations. For instance, 
the fission barriers and the nucleon(s) separation ener- 
gies are to a large extent related to ground-state nu- 
clear deformations. Theoretically, these low-lying nuclear 
structure phenomena can be well described by the mean 
field approach. Though the self-consistent Hartree-Fock 
or relativistic mean-field microscopic theories have been 
well developed during the past several decades, the well- 
known Nilsson or Woods-Saxon (WS) one-body mean- 
field potentials are still widely used in nuclear structure 
research and never stop developing. We have conducted 
some related studies [30-36] by using the PES calcula- 
tion method within the framework of the macroscopic- 
microscopic models. 


The present investigation focusing on the ground-state 
properties (e.g., shapes and binding energies) of even- 
even nuclei along the valley of stability up to the pre- 
dicted superheavy island (near the ?9§Flig4 nucleus) are 
performed in term of the realistic WS potential (a good 
approximation of the self-consistent Hartree-Fock mean- 
field), combining two different macroscopic and micro- 
scopic parts. Such an empirical potential poses a high 
predictive power while its computing algorithms are rel- 
atively simple. Very recently, we performed the investiga- 
tion of the propagating uncertainties from the WS poten- 
tial parameters [36]. However, the single-particle mean 
field is not of our primary concern here, from the point of 
view of model the main contribution of the present work 
is to provide a set of new macroscopic-model parame- 
ters by the least-square fitting to recently experimental 
data [10]. In addition, from the special footholds, we 
exhibit and analyze the experimental half-lives and two- 
nucleon separation energies of available even-even nuclei. 
With the improved PES method, we evaluate the even- 
even nuclei up to the predicted “island of superheavy ele- 
ments” along the G—stability line. A crudely bound limit 
based on two-nucleon separation energies is proposed by 
an empirical extrapolation which does not include any 
interaction assumption. 

The paper is organized as follows. In Sec. I, we 
briefly introduce the unified procedure of two potential 
energy surface methods used and simultaneously provide 
the necessary references. The results and discussions are 
presented in Sec. III. At last, we give a summary in Sec. 


IV. 


II. THEORETICAL FRAMEWORK 
A. General description of PES method 


For a long time many results have been obtained by the 
PES method with the realistic WS potential. A short 
presentation of the unified procedure of this treatment 
will be given below, including the necessary references. 
The basic idea in this approach is that the total binding 
energy of a nucleus can usually be decomposed in two 
parts, 


Erota (Z, N, 8) = Emac(Z, N, 8) + Emic(Z, N, 8) (1) 


where Emac is the macroscopic bulk-energy term, being 
a smooth function of Z, N and/or deformation, and 
Emic is the microscopic quantum corrections calculated, 
e.g., from a phenomenological single-particle potential 
well. 6 denotes a set of deformation coordinates in the 
deformation space (82, y, (4). In general, the calculated 
procedure is carried out in five steps, e.g., as mentioned 
in Ref. [37]: 


Specify nuclear shape (po- 
tential) determined by a set 
of deformation parameters 


| 


Generate the single-particle po- 
tential for protons or neutrons at 
each selected deformation point 


| 


Solve the static Schrodinger equation to 
obtain the eigenvalues (single-particle lev- 
els) and eigenfunctions (wave functions) 


J 


Calculate microscopic (shell and 
pairing) quantum corrections 


| 


[Calculate macroscopic (e.g.,liquid-drop) energy 


The total potential energy can obviously be given by the 
sum of macroscopic and microscopic energies obtained in 
the last two steps. Some different variants of the PES 
method have been developed according to the different 
ways of treating these five steps. 

In this work, we describe the nuclear shape by means of 
the parametrization of the nuclear surface with the help 
of spherical harmonics. The nuclear potential is calcu- 
lated in a multi-demensional space (62,7, 84) including 


the axial and non-axial deformation degrees of freedom. 
The single-particle energies (eigenvalues) and the corre- 
sponding wave functions (eigenfunctions) are obtained 
by solving the static Schrodinger equation with a WS- 
type nuclear potential, including a central potential, a 
spin-orbit coupling potential and the Coulomb potential 
for the protons (e.g., see Sec. 2 in Ref. [38] for detailed 
expressions). That is, by giving a set of the WS param- 
eters (e.g., the Universal parameter set is adopted here; 
see Table 1 in Ref. [38]), the deformed WS potential can 
be generated numerically at each (82, y, 84) deformation 
lattice. The WS Hamiltonian matrix elements are calcu- 
lated by using the basis of the axially deformed harmonic 
oscillator in the cylindrical coordinate system. Then the 
single-particle solutions (e.g., levels and wave functions) 
are given by diagonalizing the Hamiltonian matrix. Next, 
by using the single-particle levels obtained, the shell and 
pairing corrections at each deformation point (62, y, B4) 
are respectively calculated in terms of the Strutinsky 
method [39] and Lipkin-Nogami (LN) method [40]. It 
is worth noting that the LN method allows avoiding the 
spurious pairing phase transition which may occur in the 
calculation of the traditional Bardeen-Cooper-Schrieffer 
(BCS) type. The monopole pairing is taken into ac- 
count and its strength, G, is estimated by the average 
gap method [41]. The adopted pairing windows for both 
protons and neutrons contain dozens of single-particle 
levels, e.g., half of the particle number Z and N, below 
and above the Fermi surface. 

The macroscopic energy is usually calculated by the 
liquid drop (LD) model. Nowadays, there exist many 
phenomenological LD models in the literatures. Below, 
we list several frequently used ones in the macroscopic- 
miscroscopic calculations, such as the standard LD 
model [42], the finite-range droplet model [43], the 
Lublin-Strasbourg drop model [44], etc. These macro- 
scopic models with slightly different properties are usu- 
ally used for calculating the smoothly varying part of the 
nuclear energy, in which the dominating terms are asso- 
ciated with the volume energy, the surface energy, the 
symmetry energy and the Coulomb energy. 


B. Two variants of PES calculations 


In the present work, we perform the calculation for 
the total binding energy (defined as negative) by two 
PES methods with slightly different treatments in macro- 
scopic and microscopic energies. In the first method, the 
macroscopic energy Ep: is calculated by a standard LD 
model [42], namely, 


N-Z 
Eni (N, Z, p) = a [i-s A 4 
taji- r]es (2) 
Z? Z? 
Bo F C4 , 


“3 4173 A 


where the model parameters a,;=-15.677 MeV, a2=18.56 
MeV, c3=0.717 MeV, cy=-1.2113 MeV, and «=1.79. The 
relative surface and Coulomb energies B, and Be are only 
functions of nuclear shape. The shell and pairing cor- 
rections are considered in the microscopic energy Emici 
which reads [39], 

Emici(Z, N, B) = Esneu(Z, N, b) + Epair(Z, N, b). (3) 
Note that the treatments for the microscopic corrections 
can be found in Refs. [39, 40] (also, cf. Refs. [30-35] and 
references therein). Since these are widely used and rel- 
atively standard, it is not necessary to give the detailed 
description here. As described above, the total binding 
energy (in order to compare with data, the absolute value 
is used here) of a nucleus can be calculated by summing 
the macroscopic and microscopic energies at each defor- 
mation grid. The smoothed PES can finally be obtained 
by interpolating, e.g., using a spline function, between 
the calculated sampling points. The nuclear properties 
including the equilibrium deformations, binding energies, 
and so on can be analyzed already with the help of the 
minima of the PES. For short, we call this method PES-1 
in this work. 

Now let us introduce the second PES calculation which 
is similarly referred to as PES-2 below. In this method, 
the corresponding macroscopic energy is given by [45, 46], 


T(T +1 
Expa(N, Z) = av A + a,A2/3 + S T 12 


(4) 
aasymT (T + 1) Zz Zz? 

44/3 +C + Caz: 
Note that the present coefficients of Eq.(4) are obtained 
by us through a least-square fit to the binding energies of 
available 578 even-even nuclei [10]. During the parameter 
fit, the x? function is defined as 


578 
=e a (5) 


i=l 


where Bt©° stands for the calculated binding energy 
(namely, the negative Etotaų as mentioned above) and 
B:*? for experimental data. 

Table I shows the old [46] and new (present) parame- 
ters for the LD energy defined in Eq.(4), together with 
their root-mean-square (rms) deviation values (charac- 
terizing to some extent the quality of the fit). It should 
be pointed out that the rms value given by PES-1 is 
4.505 MeV, which is greater than those values obtained 
by PES-2 (e.g., 0.667 or 0.614 MeV, as seen in Table I). 
At this time, the microscopic part is defined as (same to 
that in Ref. [46]) 


Emic2(Z, N, £) =Faet(Z, N, B) (6) 
+ Esnen(Z, N, B) + Epair(Z, N, B), 


where the shell and pairing corrections Eshe and Epair 
are same to those in Eq.(3). The deformation correction 


TABLE I. Two sets of parameters of the macroscopic LD 
energy Ey pe2 [45, 46], as defined in Eq.(4). Note that the old 
parameters are taken from [46] and the new ones are refitted 
to 578 even-even nuclei in the AME2016 mass table. Two rms 
values are calculated for these 578 even-even nuclei based on 
the old and new parameters, respectively. Energies are in 
units of MeV. 


ay as asym Qasym C Ca Orms 


Old -15.707 18.302 117.481 -161.323 0.717 -0.882 0.667 
New -15.725 18.403 118.228 -165.496 0.719 -0.919 0.614 


energy is written as [38, 42, 47] 
Eae(Z, N, 8) = {[Bs(3)-1]+2x]Bo(8)—1] E? (7) 


where the spherical surface energy EO) and the fissility 
parameter y are Z- and N-dependent, cf. Refs. [42, 47]. 
The quantities Bg and Bo are same to those in Eq.(2) 
which includes the deformation energy. The summation 
of two terms E,p2 and Emic2 gives the total binding en- 
ergy in the PES-2 calculation. 

Note that the PES-1 and PES-2 calculations will give 
the same equilibrium deformations but different binding 
energies since the treatments of the deformation energy, 
shell and pairing corrections are same but the spherical 
LD energies different. With the aid of the PES calcu- 
lations, the present work will investigate some nuclear 
properties associated with the shape and binding energies 
along the 6-stability line. In addition, Xu et al [48] per- 
formed the configuration-constrained PES calculation for 
the high-K isomers where the binding energy calculated 
by the PES-1 method was used for adjusting the pair- 
ing strength based on the five-point mass formula [41]. 
Such an adjustment is of importance for reproducing the 
excitation energies of high-K states. Relative to PES- 
1 calculation, the present improvement reduced the rms 
value of the binding energies for 578 even-even nuclei 
from 4.505 to 0.614 MeV. It will be reasonable to believe 
that the improved PES calculation will give a better de- 
scription for the adjustment of the pairing strength since 
it can provide the relatively accurate binding energy (but 
this is beyond our present project and we just focus on 
even-even nuclei here). 


II. RESULTS AND DISCUSSION 


In order to evaluate the nuclear properties, some phe- 
nomenological or empirical or systematic laws are of- 
ten used since they are usually simple but valid and 
weakly model-dependent, even model-independent. As 
is pointed out by Casten et al. [49], the P (P = 


N,Nn/(Np + Nn) factor is one of the phenomenological 
quantities where Np and N, are respectively the num- 
bers of valence protons and neutrons (for more details, 
e.g., cf. Refs. [49-53]). Nuclear collectivity and defor- 
mation, especially in the heavy nuclei region, are ex- 
pected to depend on it sensitively. As is known, nu- 
clear deformation depends to an extent on the compe- 
tition between the like-nucleon pairing interaction and 
the neutron-proton interaction [49]. Such competition 
can be suitably described by the P factor defined above 
since it can be viewed as the average number of inter- 
actions of each valence nucleon with those of the other 
type. Much empirical evidence [52] indicates that n- 
p residual interaction and the pairing gap are typically 
on the order of 250 keV and 1.0-1.5 MeV, respectively. 
When each valence nucleon interacts with about 4-5 nu- 
cleons of the other type (corresponding to P ~ 4 — 5), 
the transition to deformation will occur near this do- 
main. At that time, the residual n-p and pairing inter- 
actions become competitive. Using the contour maps of 
the P factor as the background, Figure 1 illustrates the 
available stable nuclei and the 6-stability lines plotted 
in terms of the formulas N — Z = 0.4A?/(A + 200) [3] 
and Z = A/(1.98 + 0.0155A?/*) [4], together with the 
predicted drip-lines [18] and the boundaries of even-even 
nuclei synthesized experimentally. One can evaluate the 
P-factor distribution along the stability line. We can 
notice that the stability lines almost go across the top 
of the highest P-factor ‘hill’. Near this ‘peak’, the rota- 
tional bands are identified in 2°4No [54] and 7°°Rf [55, 56], 
indicating the expected collectivity and deformed prop- 
erties. 


In this paper, we choose to investigate the nuclear 
properties of 59 even-even nuclei which lie on or are the 
closest to the 6-stability line (the Green’s expression [3], 
N — Z = 0.4A?/(A + 200), is arbitrarily adopted) for 
given even proton numbers from Z = 2 to 118. Of course 
there will be 88 even-even nuclei on the stability line 
for given even neutron numbers from N=2 to 176 (the 
largest N in the observed even-even nuclei), but the ob- 
tained conclusions will be similar. As shown in Fig. 2(a), 
the plot of the P-factor as a function of Z indicates there 
are two peaks which are higher than the critical shadow 
region. To cross-check this large collective behaviour, we 
present the Ryo (= E(4})/E(27)) ratio of the selected 
55 even-even nuclei in Fig. 2(b), together with several 
critical points (see Refs. [58-66] for more details). Such 
a parameter is arguably one of the the best indicators of 
changes in low-lying nuclear structure. Though it can- 
not describe the magnitude of the deformation, it is often 
able to provide nuclear deformed properties [67]. More- 
over, these critical R4/2 values are widely used to identify 
collective properties of even-even nuclei. Figure 2(b) il- 
lustrates most of nuclei on the stability line exhibit the 
collective behaviour with Ryjz > 1.82 (the Mallmann 
critical point [58]), in good agreement with the results 
given by the P-factor plot. 


Figure 3 shows the ground-state equilibrium deforma- 
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FIG. 1. (Color online) Valley of stability and several typical 
boundaries on Z — N plane, using the P-factor contour map 
as the background. The black square symbol corresponds to a 
stable nucleus which is arbitrarily defined as one with a half- 
life longer than ~ 10° years. In such a classification, there will 
be 288 nuclides in total [1]. The gray and white solid lines rep- 
resent the two stability lines calculated by using the formula 
given in Ref. [4] and Ref. [3], respectively. The region between 
the two pink solid lines denote the nuclei have already been 
observed or synthesized [25]. The two black long-dash lines 
above and below the valley of stability respectively indicate 
the proton and neutron drip lines given by Moller et al. [57]. 
The wine short-dash line is the N = Z line, just using to 
guide the eye. 
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quantities R4/2 ratio [25] (a), P factor (b) for even-even nuclei 
along the £-stability line. 


tions 62, y and (4 obtained from the calculated PES min- 
ima for even-even nuclei along the stability line, together 
with the calculations by the folded-Yukawa (FY) single- 
particle potential and the finite-range droplet model 
(FRDM) [57] and the available experimental data [68, 
69]. As is known, the systematically collected data con- 
cerning the most important quadrupole deformation pa- 
rameter 32 can be obtained not only from reduced tran- 
sition probabilities B(2) [68, 69] but also from exper- 
imental quadrupole moment values [69]. The 82 value 
given by using the reduced transition probabilities B( £2) 
is always positive and by defintion automatically include 
the dynamic effect. The quadrupole deformation 82 cal- 
culated from the quadrupole moment values are obtained 
not considering nuclear surface vibrations since the rela- 
tionship is obtained by assuming near spherical nucleus 
shape vibrations amplitude is small in comparison to nu- 
cleus equilibrium deformation. 

Note that based on the measured static electric 
quadrupole moment Q values in the laboratory sys- 
tem [69], the 82 is calculated by [70] 


15Q  (1+1)(2I1 +3) 5 
~ 19Z(r)2 3K2 — I(I +1)’ (8) 


b2 


where J is the spin of the state and K is the angular mo- 
mentum projection on symmetry axis. Considering the 
surface diffusion properties in light nuclei (e.g., for A < 
100), the value of (r?) was evaluated using the expres- 
sion [71] 


_ 0.683 (1 + 10/3(ra0/Ro)?) 


m 1 + (mag/Ro)? i 


(9) 


where the parameters Ro = 1.07 A13 fm and ap = 0.55 fm 
(which can be obtained from the data on fast electrons 
scattering). For A > 100, (r?) = 0.6 x (1.2A1/9)?fm?. 
Usually, there are rather obvious differences between 
these two experimental 62 values extracted from mea- 
sured B(E2) and Q, especially in the soft nuclei. As we 
can see from Fig. 3(a), it seems that the difference of Bo 
values from two different kinds of measurements supports 
to consider the vibrating effects in the light nuclei during 
the calculations. In the heavier nuclei, our calculations 
agree with the available data and the results given by 
Möller et al [57]. Moreover, the calculated deformations, 
e.g., G2, are in good agreement with the emperical anal- 
ysis, e.g., cf. Fig. 2. The y and 84 deformations are at- 
tracting much interest in the research of nuclear shape as 
well. The present calculations will provide some useful 
information but the anomalous points (e.g.,the anoma- 
lous 62 in 12°Dyo4, as seen in Fig. 3(a)) will need to 
verify. 

Nuclear stability (or rate of decay), which is governed 
by a combination of quantum mechanical rules, nuclear 
forces, and electrostatic charge, is usually measured in 
terms of half-life. In Fig. 4, we present the half-lives of 
59 isotopic chains, from Z = 2 to 118, taking the selected 
nuclei on the (6-stability line as footholds. That is, the 
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(Color online) Calculated ground-state equilibrium deformations (a) 82, (b) y and (c) 64 for even-even nuclei with 


10 < Z < 118 along the -stability line, compared with the FY+FRDM (FF) calculations [57] and the available data [68, 69]. 
Note that the selected nucleus (e.g., N or A) for a given Z is determined by the Green’s expression [3], N— Z = 0.44?/(A+200), 
and two PES calculations give the same equilibrium deformations. See the text for more details. 


half-lives are plotted against N but the corresponding N 
of nuclei on the stability line are used for the reference 
coordinates. After such coordinate shifts, one can easily 
find that the half-lives in each isotopic chain systemat- 
ically decrease as the neutron number N moves away 
from the reference point at least up to the Pb isotopes. 
Above Z = 82, it seems that the half-life properties loose 
the systematic behaviour visible for the lighter nuclei. 
From the present data with Z > 82, regardless of the N 
number, no stable nuclei exist except for several so-called 
primordial radioactive nuclides (e.g., 732Th, ?35238U and 
244Pu), 

Moreover, one can notice that two maxima may ap- 
pear during the half-life evolution though both are un- 
stable. In addition, it is found that below the Z = 50 
Sn isotopes there are more data on the right side of the 
stability line while after this isotopes more data are ob- 
served on the neutron-deficient side. In particular, all the 
observed data locate on the left side of the stability line, 
namely the neutron-deficient direction. To explore along 
the stability line, the neutron-rich projectile-target com- 
binations seem necessary. Note that the ?940Os176 nucleus 
with the largest Z up to now has been identified but there 


is neither half-life information available nor two-nucleon 
separation energies. It is certainly expected that the half- 
lives will rapidly increase as approaching the 77$Flig4 nu- 
cleus, verifying the prediction of the ‘island of stability’. 
Anyway, from the new angle of view (standing on the 
G-stability line), it seems that one can easily notice some 
systematic properties which may be helpful for deducing 
the related physics behind them. 


The binding energy and the nucleon(s) separation en- 
ergy are of importance for nuclear structure. They can 
provide important information on the relative stability of 
nuclei, the nucleon-nucleon interactions in the nucleus, 
the bound limits of nuclei, the shell gaps, and so on. [72]. 
Under the same horizontal-coordinate convention to Fig- 
ure. 4, Figure 5 illustrates the available two-nucleon sep- 
aration energies for 59 even-even isotopic chains. As is 
well known, separation energies of isotopic (isotonic) nu- 
clei of a given parity type (e.g., even-even, even-odd, 
odd-even, or odd-odd) follow linear systematics within 
each shell region if plotted in a function of N (Z) [10]. 
Indeed, as seen in Fig. 5, one can see such linear sys- 
tematics: the two-neutron separation energy generally 
decreases with increasing N, reflecting the symmetry en- 
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FIG. 4. (Color online) The experimental half-lives (in log,)-scale) of even-even nuclei along 59 different isotopic chains, for 
each subfigure, taking the corresponding neutron number N (cf. the dash line and the N number given for each isotopic chain) 
of the even-even nucleus on the -stability line as a reference. The data are taken from the Ref. [25] and the red solid circle 
denotes the lower limit of data (e.g., for 140 Te Tı/2 > 300 ns, the data 300 ns is adopted). 
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FIG. 5. (Color online) Similar to Fig. 4, but for two-proton and two-neutron separation energies. The data are taken from 
the Ref. [10]. 


ergy and the shell filling, while at certain values of N the 
slope of the curve changes sharply [72], corresponding to 
the shell closure; the opposite trend for two-proton sep- 
aration energies. Meanwhile, it can be found from Fig. 5 
that two lines for two-proton and two-neutron separation 
energies in each isotopic chain cross at the stability line 
and the absolute values of the slopes of these two lines 
have a decreasing trend with the increasing proton num- 
ber Z. Above 763Noi62, experimental data are scarce but 
the evolution trends seem to be valid and await experi- 
mental clarification. 
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FIG. 6. (Color online) Comparison of discrepancies be- 
tween experimental and calculated (a) binding energies and 
(b) binding energy per nucleon based on the first (PES-1) and 
the second (PES-2) methods (e.g., see the text for more de- 
tails) for even-even nuclei along the -stability line. Proton 
number 12 is the lowest value of Z that we consider in the 
present investigation. Experimental data are from Ref. [10]. 


Figure 6 shows the comparison of discrepancies be- 
tween experimental and calculated binding energy and 
binding energy per nucleon based on two PES methods 
for even-even nuclei along the {-stability line. The dis- 
crepancies of binding energy and binding energy per nu- 
cleon are defined as AB = B&?-(Z, N) — Btbeo-(Z, N) 
and AB/A = B°?-(Z,N)/A — Bthe°-(Z, N)/A, respec- 
tively. Obviously, one can notice that the improved 
method gives better agreement with experimental data, 
especially for nuclei above N = 40. The improvement in 
the second method reduced the root mean square error in 
the calculated binding energies from 4.305 MeV to 0.533 
MeV and the error in the calculated binding energies per 
nucleon from 0.030 MeV to 0.010 MeV. 


2p 


AS, (MeV) 
UW ER © 
=] Ur = 


j= 
Nn 


0 20 40 60 80 100 


FIG. 7. (Color online) Similar to Fig. 6, but for (a) two- 
proton seperation energies Sz, and (b) two-neutron seperation 
energies S2n. Data are taken from Ref. [10]. 


Figure 7 shows the discrepancies between experimen- 
tal and calculated two-nucleon separation energies, e.g., 
defined by ASen = SP (Z, N) — Sth°-(Z, N) for neu- 
trons and S2,(Z,N) = B(Z, N) — B(Z, N — 2) (similarly, 
for protons). It should be pointed out that we take the 
convention that binding energy B(Z, N) is positive. As 
seen in Fig.7(a), the differences AS, have a complicated 
zigzag pattern in these two calculated results. This be- 
haviour may be due to the presence of large shell effects 
at proton closures which lead to jump in the binding 
energy. However, it is interesting that such behaviour 
seems not to occur for neutrons, as seen in Fig. 7(b). 
One can find that the discrepancies between theory and 
experiment have a decreasing trend with increasing Z. 
The improvement in the second method reduced the root 
mean square error in the calculated two-proton separa- 
tion energies from 0.966 MeV to 0.926 MeV and the error 
in the calculated two-neutron separation energies from 
0.802 MeV to 0.368 MeV. 

Based on the improved PES method, Figure 8 shows 
the calculated two-proton and two-neutron separation 
energies for even-even nuclei up to the predicted super- 
heavy island along the stability line and their differences, 
together with the comparison with experiments. It seems 
that the two-neutron separation energies against Z rep- 
resent an approximately linear relationship, decreasing 
with increasing proton numbers. It is known that two- 
neutron separation energies approximately equal twice 
the chemical potential A. As illustrated in Fig. 8(a), 
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FIG. 8. (Color online) (a) Calculated two-nucleon sepa- 
ration energies S2, (circles) and Sən (squares) of even-even 
nuclei along the stability line against Z up to 118 in terms 
of the second PES method; the results of the linear least- 
square fit are also shown, where r is the Pearson correlation 
coefficient. (b) Comparison of differences between experimen- 
tal (solid squares) and calculated (using the improved PES 
method, cf. open circle symbols) values of S2n-S2p. Data are 
taken from Ref. [10]. 


the two-nucleon separation energies are still more than 
8 MeV for the nuclei up to Z = 118 on the stability line, 
indicating that the nucleons will be strongly bound at 
these zones. Based on a simple linear least-square fit, we 
can crudely estimate the bound limit (e.g., when S = 0, 
Z = 202 +8) along the stability line. Combining the for- 
mula of the -stability line, it is easy to determine that 
the heaviest even-even bound nucleus may be (or be close 
to) the 375202 nucleus (between 344194 and 98210 along 
the stability line). From Fig. 8(b), one can see that both 
theoretically and experimentally, the differences between 
Son and Səp exhibit an obvious staggering behavior but 
the calculated results agree with available data very well, 
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IV. SUMMARY 


In summary, we have systematically investigated the 
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line based on the PES method in the (62, y, 84) de- 
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ergy) conditions in the future. 
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